Self-consistent calculations of loss compensated fishnet metamaterials 
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We present a computational approach, allowing for a self-consistent treatment of three-dimensional (3D) 
fishnet metamaterial coupled to a gain material incorporated into the nanostructure. We show numerically that 
one can compensate the losses by incorporating gain material inside the fishnet structure. The pump rate needed 
to compensate the loss is much smaller than the bulk gain and the figure of merit (FOM = |Re(n)/Im(n) |) 
increases dramatically with the pump rate. Transmission, reflection, and absorption data, as well as the retrieved 
effective parameters, are presented for the fishnet structure with and without gain material. Kramers-Kronig 
relations of the effective parameters are in excellent agreement with the retrieved results with gain. 
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The field of metamaterials has seen spectacular experimen- 
tal progress in recent years. 1-6 Yet, losses are orders of mag- 
nitude too large for the envisioned applications, such as, e.g., 
perfect lenses, 7 and invisibility cloaking. 8 For some applica- 
tions, such as perfect absorbers, 9 the loss is not a problem. 
There are two basic approaches to address the loss. One ap- 
proach to reduce the metamaterial losses to some extent is by 
geometric tailoring of the metamaterial designs^— An effi- 
cient method to geometrically reduce the losses in metamate- 
rials is by increasing the inductance, L, to the capacitance, C, 
ratior2 and avoid corners and sharp edges in metamaterials. 13 
Another method to reduce losses is to move the real part of the 
negative index of refraction, n, away from the maximum of 
Im(n) (close to the resonance) 11 ' 12 by strongly coupled meta- 
materials. The second approach is to incorporate gain ma- 
terial into metamaterial designs. One important issue is not 
to assume the metamaterial layer and the gain medium layer 
are independent from one another^— So, there is a need for 
self-consistent calculations 2 ^— for incorporating gain mate- 
rials into realistic metamaterials. The need for self-consistent 
calculations stems from the fact that increasing the gain in the 
metamaterial, the metamaterial properties change, which, in 
turn, changes the coupling to the gain medium until a steady 
state is reached. 

A time-domain self-consistent calculation of gain in two- 
dimensional (2D) magnetic metamaterial has been recently 
reported! 2 ^! It has been demonstrated 23 that the losses of the 
magnetic susceptibility, fi, of the split-ring resonator (SRR) 
can be compensated by the gain material. So, these self- 
consistent calculations 2 ^— will guide the experimentalists to 
seek new 3D metamaterial designs where the gain medium 
will give an effective gain much higher than its bulk coun- 
terpart and reduce the losses. This large value of gain is due 
to the strong local-field enhancement inside the metamaterial 
design. With regard to experiments to reduce losses in meta- 
materials, one needs to use semiconductor gain (quantum dots 
or wells) and not use dye molecules, 25 which photo-bleach 
rapidly. Semiconductor gain enables long-term use and can 
be conceptually pumped by electrical injection. This is cru- 
cial, as applications based on optically-pumped structures do 
not appear to be realistic in the long run. However, to check 



if losses in metamaterials can be reduced experimentally, one 
can try exploratory experiments under conditions of optical 
pumping. 



In this paper, we apply a detailed 3D computational model 
to study the optical response of the fishnet structure with 
gain medium embedded in the structure. We find that com- 
plete loss-compensation is possible with gain medium and the 
FOM increases dramatically with the pump rate. 26 A similar 
time-domain calculation of gain in fishnet metamaterials has 
been recently reported. 27 The gain is described by a generic 
four-level atomic system, which tracks fields and occupation 
numbers at each point in space, taking into account energy 
exchange between atoms and fields, electronic pumping, and 
non-radiative decays. 28 An external mechanism pumps elec- 
trons from the ground state level, No, to the third level, N%, at 
a certain pump rate, r pump , proportional to the optical pump- 
ing intensity in an experiment. After a short lifetime, T32, elec- 
trons transfer non-radiatively into the metastable second level, 
7V2. The second level (N2) and the first level (N{) are called 
the upper and lower lasing levels. Electrons can be transferred 
from the upper to the lower lasing level by spontaneous and 
stimulated emission. At last, electrons transfer quickly and 
non-radiatively from the first level (A^i) to the ground state 
level (No). The lifetimes and energies of the upper and lower 
lasing levels are T21, E2 and tiq, E\, respectively. The cen- 
ter frequency of the radiation is uj a = (E2 — Ei)/h, chosen 
to equal 2tt x 1.5 x 10 14 Hz. The parameters, T32, T21, and 
T10, are chosen 5 x 1CP 14 , 5 x 10~ 12 , and 5 x 10~ 14 s, re- 
spectively. The total electron density, No(t = 0) = No(t) + 
Ni(t)+N 2 (t) + N 3 (t) = 5.0 x 10 23 /m 3 , and the pump rate, 
Tpump, is an external parameter. These gain parameters are 
chosen to overlap with the resonance of the fishnet. The time- 
dependent Maxwell equations are given by V x E = — <9B jdt 
and V x H = ee a dE/dt + dP/dt, where B = /i^ H and P 
is the dispersive electric polarization density from which the 
amplification and gain can be obtained. Following the single 
electron case, we can show 28 the polarization density P(r, t) 
in the presence of an electric field obeys locally the following 
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FIG. 1 : (Color online) Schematic of the unit cell of the fishnet struc- 
ture with the parameters marked on it. The geometric parameters are 
a x = a y — 860 nm, a z = 200 nm, w x — 565 nm, w y — 265 nm, 



50 nm, t = 30 nm and h„ 



20 nm. The thicknesses of the 



metal (silver) and gain layer are t and h g , respectively, and the di- 
electric constant of the spacer, MgF2, is 1.9. These parameters were 
used on simulations 12 and experiments.— 



equation of motion, 
d 2 P{t) , „ dP(t) 



dt 2 



r a - 



dt 



-iu 2 a P(t) = -a a AN(t)E(t), (1) 



where T a is the linewidth of the atomic transition ui a and 
is equal to 2vr x 20 x 10 12 Hz. The factor, AN(r, t) = 
N 2 (r, t) — N\{r, t), is the population inversion that drives the 
polarization, and a a is the coupling strength of P to the ex- 
ternal electric field and its value is taken to be 10 -4 C 2 /kg. 
It follows 28 from Eq. Q] that the amplification line shape is 
Lorentzian and homogeneously broadened. The occupation 
numbers at each spatial point vary according to 
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where — !— E • ^E- is the induced radiation rate or excitation 

nuj a at 

rate depending on its sign. 

To solve the behavior of the active materials in the electro- 
magnetic fields numerically, the finite-difference time-domain 
(FDTD) technique is utilized. 29 In the FDTD calculations, 
the discrete time and space steps are chosen to be At = 
8.0 x 10~ 18 s and Ax = 5.0 x 10~ 9 m. The initial condi- 
tion is that all the electrons are in the ground state, so there is 
no field, no polarization, and no spontaneous emission. Then, 
the electrons are pumped from No to N3 (then relaxing to 
N 2 ) with a constant pump rate, r pump . The system begins to 
evolve according to the system of equations above. 

In Fig. 1, we show the unit cell of the fishnet structure. The 
size of the unit cell along the propagation direction is a z . a z 
is larger than the sum of the thickness of the metallic and the 
dielectric layers 2i + s, where t and s are thicknesses of the 



metal and the dielectric layers, respectively. Notice the prop- 
agation direction is perpendicular to the plane of the fishnet 
with the electric and magnetic fields along the x and y direc- 
tions, respectively. All retrieved effective parameters are for 
this particular incident direction and field polarization. Con- 
sider two configurations, one without gain and one with gain 
(see inset of Fig. 1). In the configuration with gain, we have 
introduced two thin dielectric layers of thickness, (s — h g )/2, 
close to the metallic structure, so the gain medium will not 
be close to the metal, then quenching will be avoided. The 
dimensions of the fishnet structure 11 ' 12 are chosen such that 
the magnetic resonance wavelength at A = 2000 nm, which 
can overlap with the peak of the emission of the gain material. 
The full width at half maximum (FWHM) of the gain mate- 
rial is 20 THz and the pump rate, r pump , changes from to 
6.9 xl0 8 s _1 . 



In Fig. 2, we plot the transmission, T 



tion, R 



, (2b) and absorption, A = 1 



\t\ 2 , (2a), reflec- 
T-R, (2c) ver- 



sus wavelength for different pump rates (t and r are the trans- 
mission and reflection amplitudes, respectively). Notice the 
wavelength dependence of T and R for different pump rates 
away from the resonance wavelength, A = 2000 nm, are the 
same. Below the resonance wavelength, T increases with the 
pump rate and above the resonance wavelength, T decreases 
with the pump rate. The reflection R, below the resonance 
wavelength, decreases with the pump rate and above the res- 
onance wavelength, it increases with the pump rate. Notice 
in Fig. 2a, T without gain has a very weak resonance, and 
once the pump rate increases, the transmission clearly shows 
the resonance behavior. The same can be seen in the exper- 
iments that can use gain materials to compensate the losses 
in fishnet metamaterials. In Fig. 2c, we plot the absorption, 
A, as a function of wavelength for different pump rates. No- 
tice, as we increase the pump rate, the absorption decreases 
and finally at the pump rate of 6.6 x 10 9 s , the gain over- 
compensates the losses and the absorption becomes negative. 
In Fig. 3a, we plot the retrieved results of the real and imag- 
inary parts of the magnetic permeability, //, with and without 
gain, for normal incidence and the particular field polariza- 
tion in Fig. 1. As gain increases, the Re(/i) becomes steeper 
at the resonance wavelength and the Im(/i) becomes much 
narrower when increasing the pump rate and the losses are 
compensated by the gain material. In Fig. 3b, we plot the re- 
trieved results for the effective index of refraction n, with and 
without gain. The Re(n) becomes more negative after gain 
is introduced and the Im(n) also drops significantly close to 
the resonance. At A = 1976 nm, the Re(n) changes from - 
2.25 to -2.82 with a pump rate of 5.0 x 10 8 s _1 and the Im(n) 
drops from 1.58 to 0.54 (Fig. 3b). In Fig. 3c, we plot the 
FOM versus the wavelength for different pump rates. No- 
tice, the FOM becomes very large (of order of 10 2 ) with 
the pump rates. Comparing Im(n) slightly below the reso- 
nance at A = 1976 nm, we find the effective extinction co- 
flm(n) ~ 5.0 x 10 4 cm -1 without gain and 



efficient a 
a w 1.7 X 10 4 cm _i with gain (r pump 
Hence, an effective amplification of a = 



= 5.0 x 10 8 s- 1 ). 
-3.3 x 10 4 cm- 1 . 



This is much larger (of the order of 30) than the expected am- 
plification a « —1.3 x 10 3 cm -1 for the gain material at the 
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FIG. 2: (Color online) The transmission (a), reflection (b) and ab- FIG. 3: (Color online) The retrieved results for the real (solid lines) 

sorption (c) as a function of wavelength for different pump rates. and the imaginary (dashed lines) parts of (a) the effective perme- 
ability, /i, and (b) the corresponding effective index of refraction, n, 
without and with gain for different pump rates, (c) The figure-of- 

„ , merit (FOM) as a function of wavelength for different pump rates, 

given pump rate of 5.0 x 10" s . The difference can be ex- The gain bandwidth is 2 n THz. 

plained by the field enhancement in the fishnet metamaterial. 

There are theoretical debates^— if it is possible to ob- 
tain low loss metamaterials with negative refractive index, 

n. They have used the Kramers-Kronig (KK) relations and 5.0 x 10 s s _1 . The excellent agreement between the results 

we would like to verify that Kramers-Kronig relations work obtained from the standard retrieval method 33 (a special case 

with and without gain. In addition, we need to compare of anisotropic retrieval method 34 for normal incidence) and 

the numerically-retrieved effective permeability, fj,, shown in the Kramers-Kronig approach verifies KK relations work for 

Fig. 3a, with the calculation of fi, based on the Kramers- metamaterials coupled with gain too. KK relations are valid 

Kronig relations. In Fig. 4a, we plot the real and imaginary for passive materials and the deduction of KK relations is 

parts of the effective permeability, \i, and the results from based on the assumption that the response function of mate- 

the Kramers-Kronig relations for the pump rate of r pump = rials is analytical in the upper complex plane of frequency. 
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FIG. 4: (Color online) The real and imaginary parts of the retrieved effective permeability, (i, and the results from the Kramers-Kronig relations 
for pump rates (a) r pump = 5.0 x 10 8 s" 1 (below compensation) and (b) r pump = 6.9 x 10 8 s _1 (overcompensated). 



For strong active materials, there are poles in the upper plane 
and one must modify the KK relations by reversing the signs 
of KK relations. 32 In Fig. 4b, we plot the Re(/i) and Im(/i) 
obtained by the retrieval method and by the modified KK re- 
lations for the pump rate r pump = 6.9 x 10 8 s _1 . For this 
pump rate the Im(/i) becomes negative and we have overcom- 
pensated at the resonance frequency of 150 THz. As one can 
see from Fig. 4b, KK relations work well for strongly active 
materials. 

In conclusion, we have proposed and numerically solved 
a self-consistent model incorporating gain in the 3D fishnet 
dispersive metamaterial. We show numerically that one can 
compensate the losses of the fishnet metamaterial. We have 
presented results for T, R, and A without and with gain for 
different pump rates. Once the pump rate increases, both T 
and R show a resonance behavior. We have retrieved the ef- 



fective parameters for different pump rates and the losses are 
compensated with gain. Kramers-Kronig relations of the ef- 
fective parameter are in excellent agreement with the retrieved 
results with gain. The figure-of-merit (FOM) with gain in- 
creases dramatically and the pump rate needed to compensate 
the loss is much smaller than the bulk gain. This aspect is due 
to the strong local-field enhancement inside the fishnet struc- 
ture. 
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